Quantum theory of ultracold atom-ion collisions 
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We study atom-ion scattering in the ultracold regime. To this aim, an analytical model based on 
the multichannel quantum defect formalism is developed and compared to close-coupled numerical 
calculations. We investigate the occurrence of magnetic Feshbach resonances focusing on the specific 
40 Ca + + Na system. The presence of several resonances at experimentally accessible magnetic fields 
should allow the atom-ion interaction to be precisely tuned. A fully quantum-mechanical study 
of charge exchange processes shows that charge-exchange rates should remain small even in the 
presence of resonance effects. Most of our results can be cast in a system-independent form and are 
important for the realization of the charge-neutral ultracold systems. 



Advances in trapping, cooling, manipulation and read- 
out of single atoms and ions have led over recent years to 
a range of fundamental as well as applied investigations 
on the quantum properties of such systems. Nowadays, 
an increasing number of experimental groups worldwide 
are starting experiments with combined charged-neutral 
systems in various configurations [![. While the theory 
of atom-ion collisions is well established for high collision 
energies 0, Q , a theoretical description in the ultracold 
domain is still largely missing. 

This letter presents the first study of magnetic Fes- 
hbach resonances and the first fully quantum study of 
the radiative charge exchange process for ultracold atom- 
ion systems that includes effects of Feshbach and shape 
resonances. Here we consider only two-body collisions 
in free space, a necessary prelude to further studies in- 
corporating effects of ion micromotion or trap confine- 
ment. We develop a reliable yet manageable effective 
model of atom-ion collisions by applying multichannel 
quantum defect theory (MQDT) [J, [f|, Q based on the 
long range ion-induced-dipole potential that varies as r~ A 
at large ion-atom distance r 0, |]. This powerful tool 
has proven effective as a few-parameter approach for de- 
scribing scattering and bound states in electron-ion core 
[1] , electron-atom Q and neutral atom systems ■ Al- 
though the literature on the subject is rich, here we dis- 
cuss some details of MQDT illustrating how it works in 
the ultracold domain, so we can reveal the new and in- 
teresting ultracold ion-atom physics. We adapt MQDT 
to the atom-ion realm, utilizing the analytical solutions 
for the r -4 asymptotic potential and applying the 

frame transformation jT^, [T^] at short distances to re- 
duce the number of quantum defect parameters in the 
model. We verify the model predictions by comparing 
to our own numerical close-coupled calculations, taking 
40 Ca + — 23 Na [l3[ as a reference system. 

We describe the S-state atom and S'-state ion collisions 
with the close-coupled radial Schrodingcr equation 
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Here, n = mim a / '(m, + m a ) denotes the reduced mass, 
W(r) is the interaction matrix, and F(r) is the matrix 
of radial solutions. The wave function for N scatter- 
ing channels reads \&i(r) = J^jLi AjYj(r)Fij(r)/r where 
Yj(f) denotes the angular part of the solution for the 
scattering channel j, and the constant vector A is de- 
termined by the boundary conditions at r — * oo. The 
asymptotic channel states can be characterized by the 
hyperfine quantum numbers: fi,mf 1 and f2,nif 2 for ion 
and atom respectively, and by the angular-momentum 
quantum numbers I and mi of the relative motion of the 
atom and ion centers of mass. We denote the asymptotic 
channel states by I^q), where a = {f\f2 i mf 1 Tnf 2 lmi\. In 
the presence of a magnetic field B, the field-dressed chan- 
nel states \^> a (B)) arc linear combinations of the bare 
(B = 0) channel states \® a (B)) = E Q "^«a'(-B) 
In the asymptotic channel basis the interaction matrix is 
diagonal at large distances 
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E - W(r)] F(r) = 0. 
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where Ef° are the threshold energies for the channel i 
including the hyperfine energies and Zccman shifts, and 
C4 = ae 2 /2 with a denoting the static dipolar polar- 
izability of the atom and e is the ion charge. At typi- 
cal distances Rq where the short-range exchange interac- 
tion takes place, the interaction matrix becomes diagonal 
in the IS representation, characterized by total electron 
spin S = Si + S2, the total nuclear spin I = ii + 12, the 
total hyperfine angular momentum F = fi+f2 = I + S, 
and its projection Mp on the axis of quantization, where 
Si, S2 are electron spin of ion and atom respectively, and 
ii, \2 denote their nuclear spins, respectively. We label 
those channels by (3 — {ISFMplm{\. 

With the long-rang e atom-io n one can associate the 
length scale R* = \j2C^jW- and energy scale E* = 
h 2 / [2(i(R*) 2 ]. The length R* can be also related to 
the position of the last node of the zero-energy s-wave 
radial wave function for the scattering length a = 0: 
L = 2/ttR*. The maxima due to the centrifugal bar- 
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\)R* and have heights 
Hence E* determines the contri- 



rier occur at r max = 
of£ max = |Z 2 (Z+1) 2 
bution of higher partial waves in the atom-ion scattering. 

MQDT. The basic idea of MQDT ^MM is the sepa- 
ration between the short-range and the long-range prop- 
erties of the scattering wave function. The short-range 
wave function is insensitive to the total energy E, and in 
the case of atom-ion collisions also to the relative orbital 
angular momentum [hi ]. On the contrary, at large dis- 
tances the wave function vary rapidly with E, as well as 
other scattering quantities defined at r — » oo. By intro- 
ducing a set of short-range quantum-defect parameters, 
MQDT allows one to determine all scattering and bound- 
state properties in a wide range of collision energies. 

In MQDT one starts by choosing a set of reference po- 
tentials {Vj(r)}, that reproduce the asymptotic behavior 
of the interaction matrix at large distances: Vj (r) r — > 
Wjj(r), but otherwise can be arbitrary. With the refer- 
ence potentials Vi (r) one can associate a pair of linearly 
independent solutions fi(r) and §i(r) 
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fi(r) = cti(r) sin/3i(r), 
<7i(r) = oti(r) cos Pi(r). 



(3) 
(4) 



The amplitude cn(r) fulfills the inhomogeneous Milne 



equation: [d 2 /dr 2 + fc,(r) 2 ]a(r) 
local wavevector fcj(r) : 
phase dpi/dr = I /a 2 . 



a 



! (r) [15f, with the 



= y/2fi(E -Vi(R)/h, while the 
We impose WKB-like bound- 
ary conditions at small distances: oti(r) = 1 / -\Jki (r) , 
aj(r) = 0, which makes the functions fi(r) and (ji(r) 
weakly dependent on energy. The exact solution to 
Eq.([T]) at large distances (where Vj(r) = Wjj(r)) can 
be expressed in terms of fi(r) and <?;(V): F(r) r -^3 > 
'f(r)+g(r)Y(£)] A. Here, f(r) = diag^r)], g(r) = 

diag[5i(r)], and Y(E) is the energy insensitive (Y(E) = 
Y) quantum-defect matrix that plays a central role in 
the quantum-defect analysis. 

Analytical solutions. The Schrodingcr equation with 
the polarization potential C4/V 4 can be solved analyti- 
cally, in terms of Mathieu functions of imaginary argu- 
ment [ll|, [l6[ . For our choice of MQDT reference poten- 

E°°, the solutions: 



tials: Vi(r) = C4/V 4 + 1(1 + 
f, g arc singular at r — > (Ti 
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f(r) = rsin(—R*/r + tp), 
g(r) = r cos(— R* /r + tp), 
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Here, y denotes some short-range phase, that can be re- 
lated to the s-wave scattering length a = —R* cot tp of the 
C 4 /r 4 potential [13 ■ In view of the arbitrariness of the 
reference potential, we set tp = 0, which does not affect 
the physical scattering matrices, that are not sensitive to 
the short-range behavior of the reference potentials. 

Frame transformation. The frame transformation 
(FT) is a unitary transformation between channels a and 
0'- U a /3 — (fx f 2171^171 f 2 \ISFMp), that can be expressed 
in terms of the usual Clcbsch-Gordan and Wigner 9j 
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FIG. 1: (Color online) Elastic collision rate in the singlet 
channel J 4 1 S + (upper panel) and rate of radiative A 1 E + — > 
X 1 E + charge transfer (lower panel) versus collision energy. 
Different theoretical approaches are compared (see text). 
Both the elastic and the total inelastic rate are decomposed 
on selected partial waves. 



symbols (see e.g. [l(|). The application of FT requires 
i?o R*, which is typically very well fulfilled in the 
atom-ion collisions. At distances r > Rq the exchange 
interaction is no longer present, and for r <JC R* we can 
safely neglect both the centrifugal potential and hyper- 
fine splittings, and also ignore the higher order dispersion 
terms: Ce/r 6 ,Cs/r 8 that lead only to small corrections 
(see below). Then for R Q <r<$:R*, W t] (r) = S^C^/r 4 , 
and the wave function in each of the channel is simply 
given by the linear combination of the two short-range 
solutions © and ©, both in the asymptotic and (IS) 
representations . 

In view of the assumed form of the reference poten- 
tials, and the choice tp — 0, the quantum-defect matrix in 
IS representation takes the form Y^, = Sppi [as(p)]^ 1 , 
where S(j3) is the total electron spin in the channel (3, 
that takes values S((3) = 0, 1 corresponding to the singlet 
a s and triplet a* scattering lengths, respectively. Then, 

Y in the basis of the asymptotic channel states reads 

Y = Z(S)UY( /s 'U t Zt( J B) where 1(B) is the transfor- 
mation from bare to the magnetic-field-dressed states. 
Finally, the S and K matrices describing all the scat- 
tering properties including the cross-sections, are calcu- 
lated directly from Y with the help of the quantum-defect 
functions [5|, |6( , that are determined by the asymptotic 
behavior of the short-range normalized solutions f,g. 

Results. We now apply our model to a system of 40 Ca + 
and 23 Na for which the long-range parameter C4 = 81.35 
a.u. and the short-range ab-initio potential curves arc ap- 
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proximately known [13j. The close-coupled Schrddingcr 
equation (JTJ) is solved both within our MQDT approach 
and exactly using standard numerical methods. 

Unfortunately, potentials calculated by ab-initio meth- 
ods are usually not sufficiently accurate to predict the 
value of the associated scattering lengths. Therefore, for 
our calculation we will either assume typical magnitudes 
of a s and at on the order of R* or we vary them. It is 
interesting to note that for a r~ 4 potential the probabil- 
ities to find a positive or a negative scattering length are 
equal in contrast to van der Waals r~ e interaction where 
their ratio is 3:4 [H[. 

The 40 Ca ion has vanishing nuclear spin (ii = 0), 
whereas the 23 Na atom has 12 = 3/2, resulting in a to- 
tal hypcrfinc angular momentum /a = 1 and f% = 2. 
Neglecting small anisotropic spin interactions, rotational 
invariance implies that I, m;, and Mp are conserved 
quantities. In our calculations we will consider colli- 
sions in the Mp = 1/2 block, and use the FT to ob- 
tain Y in the asymptotic channel representation. At 
B = the off-diagonal elements of Y are proportional 
to l/a c = l/a s — l/a f , where a c is an effective scattering 
length characterizing the strength of channel coupling. 

We begin our analysis by considering a single 
potential as entrance channel. The upper panel of Fig. [T] 
shows the rate vai of elastic Ca + + Na collisions as a 
function a collisions energy for different partial waves, 
with v the relative atomic velocity and o~i the partial 
wave cross section. The calculation is performed us- 
ing a potential parameterized by a scattering length of 
a s = — 0.961i?* = — 2000ao. One can observe several 
dips of the collisional rates, corresponding to Ramsauer 
minima. We note that because of the relatively low cen- 
trifugal barriers, partial waves with I > give a signif- 
icant contribution already in the /xK regime. At small 
energies the rates for I > behave asymptotically as 
aiv ~ E 3 / 2 , in agreement with the threshold behavior 
typical of r~ 4 potentials [l!|. For I = 4 and at higher 
energies some discrepancies between the numerical and 
MQDT approach can be observed, resulting from break- 
down of the assumption of angular-momentum insensitiv- 
ity of the quantum-defect matrix in the analytical model. 
On the basis of the MQDT model, one can argue that the 
elastic collisional rates in the triplet channels, would be 
of the similar order, with the particular structure of max- 
ima/Ramsauer minima determined by the value of at- 

The lower panel of Fig. [1] shows the total rate Wt r of 
inelastic A 1 Y, + — > X 1 Y> + collisions associated with ra- 
diative charge transfer (see [l3[ for details). The numer- 
ical results (first three entries in the legend) have been 
obtained in the distorted-wave Born approximation by 
summing contributions from all free-free and free-bound 
transitions [2(j. Note that molecular ions are formed 
with very high probability (~ 96%) in the electron ex- 
change process. The figure also shows the total charge- 
exchange rate obtained by a computationally much sim- 
pler single channel calculation with a position-dependent 
Einstein coefficient A{r) (see e.g. 2l|). In addition we 
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FIG. 2: (Color online) Field dependent s-wave scattering 
length (upper panel) , charge transfer rate for s-wave collisions 
at E = InK (middle panel) and energies of I = bound states 
(lower panel) as a function of the magnetic field strength B, 
for a s = —0.9R* and at = —a 3 . Close-coupled numerical cal- 
culations with (solid line) and without higher-order dispersion 
terms (dotted line), and the MQDT results (dashed line) are 
compared. Results from different approaches are undistin- 
guishable on the scale of the bottom panel. 



show the result of our analytical model where the charge- 
transfer probability has been described in the semiclas- 
sical theory, while MQDT provides a proper rescaling 
of the wave function between reaction and asymptotic 
zones. One can observe that the three numerical and an- 
alytical approaches give virtually the same result for the 
total (free-free plus free-bound) charge- exchange rate at 
low energy It is only at higher energies, in particular 
near the peaks associated to shape resonances, that the 
MQDT model becomes less accurate. This can again re- 
sult from the modification of the quantum-defect matrix 
Y for higher I. 

In contrast to the purely scmiclassical result of Ref. 
[l3[, both elastic and inelastic quantum cross sections 
appear to be strongly affected by the presence of shape 
resonances up to mK regime; see Fig.[TJ We have verified 
that this is a general behavior not related to the partic- 
ular choice of the singlet scattering length assumed for 
the calculation. The influence of shape resonances can 
also be observed in the thermally averaged rate which 
only begins to saturate to the asymptotic value in the 
mK regime. 

Next, we take into account multichannel hypcrfine ef- 
fects and the presence of an external magnetic field. 
Fig. [2] shows the variation of the effective scattering 
length a versus magnetic field for sample values of the sin- 
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FIG. 3: (Color online) Contour plot of the dimensionless 
quantity (2/tt) arctan(a/i?*), versus singlet and triplet short- 
range phases 4> a — — arctanl/a s and (j>t = — arctanl/at and 
magnetic field B. In the leftmost panels 4> 3 is fixed, in the 
rightmost panels <f> t is fixed. 

glet/triplet scattering lengths: a s = — 0.9R* = — 1873ao 
and at = — a s . Several broad magnetic resonances can be 
observed in the experimentally accessible range of mag- 
netic fields. As for neutral atoms, such resonances can in 
principle be used in practice to tune the atom-ion effec- 
tive interaction. 

The MQDT results agree very well with the numerical 
results without higher-order dispersion terms, whereas 
one can observe very small discrepancies in the magni- 
tude of a in comparison to the calculations performed 
with the full potential. Each resonance is associated 
with a bound molecular level crossing the energetically 
lowest atomic threshold, as shown in the lowest panel. 
The width of a resonance depends on particular molec- 
ular state crossing the threshold. For small interchan- 
nel coupling l/a c bound states occur mainly in a single 
asymptotic channel a, whereas for larger l/a c and weak 
magnetic fields, they are better characterized in terms of 
hyperfmc angular momenta F, /i, /2- 

The middle panel compares the radiative charge-trans- 
fer rate calculated numerically with the closed-coupling 
method and a local Einstein coefficient, and MQDT with 
the semiclassical description of the short-range charge- 
transfer process. In the vicinity of resonances the charge- 
transfer rate is strongly enhanced as expected, and the 
MQDT model seems to slightly overestimate it. 



The range of magnetic fields in which the first Fcsh- 
bach resonances should occur can be estimated from the 
energy location E°° of the associated bound states for 
a = oo and B = 0. For a pure r~ 4 potential such en- 
ergies are E°° / E* = 105.8,1180,5208,... Letting A^i be 
the difference of magnetic moments between the Fesh- 
bach molecule and free atoms and assuming only linear 
Zeeman shift for small values of B, we expect the first 
resonances to be found at B < B max = E°° /Afi, with 
B max [mT] = 0.086,0.962,4.26 for Na-Ca+ system. 

In order to investigate the dependence of the position 
and strength of Feshbach resonances on the interatomic 
potentials, we show in Fig. [3] the variation of the quan- 
tity 2 arctan(a/i?*)/7r with the singlet and triplet short- 
range phases (p s> t = — arctan l/a Sj t and magnetic field B. 
Please note that typically several resonances should be 
observable below B ~ 20 mT. One can remark that the 
resonances are usually relatively broad and only become 
very narrow when a s and at are comparable, and the in- 
tcrchannel coupling 1 / a c is small. Apart from the vicin- 
ity of resonances the field dependent scattering length is 
mainly determined by the value of a t , that again is in 
turn related to the structure of the Y matrix. 

In conclusion, a relatively simple quantum-defect mo- 
del accurately describes atom-ion collisions in the ul- 
tracold domain, as verified by comparison with numer- 
ical coupled-channel calculations for the reference sys- 
tem 40 Ca + and 23 Na. Our calculations predict that sev- 
eral magnetic Feshbach resonances should be available at 
relatively small values of magnetic fields to control the 
atom-ion interaction. Radiative charge exchange rates 
remain relative small even in the presence of Feshbach 
and shape resonances so that elastic collisions can be con- 
trolled without introducing unwanted losses. 
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